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Recently, Dorey and Tateo have investigated functional relations among Stokes 
D \ multipliers for a Schrodinger equation (second order differential equation) with a 

"^ ' polynomial potential term in view of solvable models. Here we extend their studies 

to a restricted case of n + 1— th order linear differential equations. 
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1 Introduction 

In remarkable papers]!], §], Dorey and Tateo find marvelous connections between a 
Sclirodinger equation with a polynomial potential term and solvable models. The suc- 
cess integrates two ingredients. The exact WKB analysis0-[|^ is on the one side. The 
development in studies of solvable models and of the solvability structure in (perturbed) 
CFT in], ]12], ]T3], ^M via the thermodynamic Bethe ansatz|]l6], ]r7| and nonlinear integral 



equations is on the other |T^-[^ . 



In view of the eigenvalue problem in quantum mechanics, Ti"^ = E^^ the quantity of 
interest is the spectral determinant, det(7i — E). For the potential term x"^^^ with M = 
integer, it coincides with the "vacuum expectation value" of a particular member of a 
fusion transfer matrices in a certain field theory possessing Uq{sl2) [|], ^. The rest of 
members in the fusion hierarchy are identified with the Stokes multipliers and their gen- 
erahzation in the vacuum sector. Thus the result provides a unified view of Stokes 
multipliers and spectral determinant. 

The spectral determinant for the wider class of potentials x^*^ + i{i + l)/x^ with M 
general can be treated within a framework using Baxter's Q— operator and the quantum 
Wronskian relation [Q. 

We extend a part of studies in to the higher order differential equation case, 

<9"+^y+(-l)'^(a;'" + A"+^)t/ = 0, (1) 

with m integer. 

The A = case is essentially equivalent to the Turrittin equation of which solutions are 

calculated in terms of Meijier G— functions. There have been several results on this case 



11,11, HI. 



We consider the Stokes multipliers associated with eq. (]!]). They satisfy a set of func- 
tional relations in the complex A plane. We will show that nontrivial solutions to relations 
are expressible by the quantum Jacobi-Trudi formula which appears in fusion transfer 
matrices related to UqiA'n^) p8], H. 



We note the very recent result in |]5D[ where the eigenvalue problem of eq. (]!]) is argued 
for n = 2 case with x G [0, -|-cx)] though a related but a different method . 

2 Asymptotic Expansion and Stokes Coefficients 



We refer |3T]] to readers for the background on the subject. 

Let us first discuss the asymptotic behavior of a slightly generalized differential equa- 
tion, 

9"+^y + (-l)"P(x)i/ = 

m 

P{x)=J2<^,x--' (2) 

j=0 



where aj are complex numbers and oq = 1. P{x) will be referred to as the potential 
term in analogy with Schrodinger equations. [| The factor (—1)" is not essential. It can be 
adsorbed into re-definition of the angle of x and aj>i. For later convenience, we include 
this factor throughout this report. 

Now that a; = cxo is an irregular singular point of the equation, analytic properties 
of solutions depend on sectors in the complex x plane. Let Sk be a sector in the plane 
satisfying 

o 

largx — k6\ < - 
for X & Sk, where 9 = ^l^+i • ^^ then study the asymptotic behavior of a subdominant 



solution in Sq. Following |Q 0], we define bh{h = 1, 2, ■ ■ • ) by the relation. 



(1 + ^ afcX-'=)^/("+^) = 1 + 5] ^hX-^ 
k=i h=i 

A key function E{x, a) is determined by bh, 

{l + ^hhX-^)x'^'^''^^^dx 
h=l 



I 1 "^ h 



m + n + 1 ^^ ^-h + l 

h=l n+1 

where hm = N ior m = N{n + 1) — j, (j = 1, ■ ■ ■ , ra). Here a stands for (ai, a2, ■ • ■ , am) 
In addition, we introduce an exponent Um by 

_ ^, for m 7^ mod n + 1 , . 

J^m- z-n + (^ + l)bh^+i, for m = mod n + 1 . ^"^^ 

Then we have a theorem. 

Theorem 1 In Sq, there exists a subdominant solution y{x, a) to eq.l^ which admits the 
following asymptotics, 

y{x,Si) ~ 3;--™/(n+l)g-S(x,a)^ 

The range of the existence of the asymptotic expansion (|) (forgetting the subdomi- 
nance) actually extends over Sq. The precise determination of the range, however, poses 
a nontrivial problem even for n = 1 case. See elaborate discussions in sections 9 and 10 

^ In the case of singularity at infinity of rank 1, i.e., P{x) — pq + pi/x + • • ■ in the above, asymptotic 
solutions and numerical algorithms are discussed in p3| for a wider class of higher order differential 
equations. 



of IQ . The straightforward extension of the argument may lead to the conclusion that 
y{x, a) admits the asymptotic expansion (Q) in an open sector, 

n + 2 



argx < 



-IT. 



(5) 



m + n + 1 

See the appendix. 

The intriguing feature in the differential equation (§) is the following symmetry in 
rotating the x— plane. 

Theorem 2 Denote the above solution byy{x,Si). Then 

l/fc(x,a):=y(xg-^G«(a))g"'=/2 (6) 

is also a solution to eq.iji^. 

The parameter q signifies exp(z6') = exp(i ^^^ ). The operation G^'^^a) is defined by 
GW(a) = G(i)(G('=-i)(a)),A; > 2 and G^^\af={ai/q,a2/q'^,---cim/q'^)- 

A fundamental system of solutions (FSS) in S^ is formed by {yk, yk+i, ■ ■ ■ , Uk+n)- This 
is shown for the Turrittin equation in [^. For general, it is not simple to show this 
directly. We use the range of the validity of asymptotic expansion in eq.(|) and justify 
the linear independence as follows. 

We introduce a (n + 1) x (n + 1) matrix ^k{,x) and the Wronskian Wk := det^k{,x) 



^kix) 



( Vk, 
dVk, 



2/fc+i, 



Z/fc+n \ 

dyk+n 



(7) 



\(9"?/fc, <9"?/fc+i, • ■ ■ , d'^yk+nj 

Then eq.(|^) implies that yu+ji 3 = 0, ■■■ ,n has the asymptotic expansion (H) in the 
closed sector Sk+(n+i)/2 U <Sk+{n-i)/2- As Wk is constant in x, one then verifies the lin- 
ear independence of these solutions by using the asymptotic expansions (^ and (|) in 
Sk+(n+i)/2 U Sk+(n-i)/2- Scc discussious in [|^ for n = 1 and 0] for n = 2 

We stress that the argument below depends on the existence of ^k{x) as FSS, rather 
than (1). 

A matrix of Stokes multipliers S^ connects FFS of Sk and Sk+i 



The linear independence of solutions fixes 5*^ in the following form. 



(8) 



S, 



(1) 



^k ' 

^(2) 
^k ' 


1, 
0, 


0, 

1, 


0, ■■ 
0, ■■ 


-, 


J-) 
^k ' 

yk 5 


0, 
0, 


0, 
0, 


0, ■■ 
0, ■■ 


-, 1 



(9) 



By Cramer's formula, one represents rjf as 



^ I Vk+l, yk+2i ■ ■ ■ 5 Vki ■ ■ ■ 5 Hk+n+l 



fc ly. 



' 9"?/fc+i, <9"yfc+2, ■ ■ ■ , ^'^l/fc, ■ ■ ■ , d'^Vk+n+l) 



(10) 



The column vector, {yk,dyk, ■ ■ ■ ,d"yk) is inserted in the j— th column in the matrix of 
the denominator. Especially, 

rf^ = {-irW,+i/W,. (11) 

Then we restrict ourselves to the case of our interest (|1]). We find that 6h„+i = and 
thus Um = nm/2 for any m. Under the operation of G^^\ G^^\am) = a^q'^ = am?"^"^- 
As ttm = A""*"^, this means G^^\\) = \q. Thus a function in the k— th sector has an 
argument Ag^. 

3 Recursion relations 

By the commensurability of the cone angle of iS^ to 27r, the product of successive m + n + 1 



Stokes matrices must be a unit matrix ||33| modulo sign factor due to the normalization 
ind). 

Regarding as a function of X, Sp. = 5'*^^^(Ag^). Thus eq.(|T^) reads, 

^(1) (Ag™+")^(i) (Ag'^+^-i) ■ ■ • S(^) (A) = (-1)"/. (13) 

The same relation can be recapitulated in terms of a generalized Stokes matrix Sf, con- 
necting $fc and $fc+£ II, |T1, 

^k = ^k+eSi'l (14) 

Similar to the above, S^ = S^^\Xq'') as a function of A. 
Eq. (|T3p is obviously rewritten as, 

^(m+n+l)^^^ = (-1)"/. (15) 

By the above definition, we have two equivalent recursion relations 

5(^+i)(A) = S^'\\q')S^'\\) (16) 

= S^'\Xq)S^^\X). (17) 

We denote the (i, j) component of S^^\X) by SlJ{X). 
The main problems in the present report are 



(1) the expression of SlJ{X) in terms of r*^"^(A)s 

(2) restrictions imposed on r^"-'(A)s by functional relations among them. 
The first recursion relation (|T^) imposes the relations, 

SiT(^) = r^'\q'^)S^i^W + Si^l,/X) {l<k<n + l) (18) 

where we formally put 5'„_^2,j('^) = 

The second recursion relation (|I^ yields, 



+ r("+i)(A)5£^,(Ag) (19) 



and 

sZt'iW = Si']{qX) {l<i<n + l,l<j<n). (20) 

The latter equation leads to 

Slti(^) = Si\''\q^X) i>j. (21) 

For i < j, 

^r(A) = ^iUg^A), (22) 

and the rhs is given directly in terms of r*^"^ . Thus we regard only S^, [ as nontrivial 
elements. 

The eq. ([T9|) can be rewritten such that it contains only S^ { in the Ihs. Let us write 



explicitly for k = 1, 

5ir(A) = r«(A)5S(Ag) + r(^)(A)5iS-^)(Ag^) + 

••• + r("+^)(A)5i'r"^(Ag"+i). (23) 

We remark properties of S{ { (A) from linear independency of FSS. The (1,1) component 
of eq. ([T^ ) for i = m + j,l < j < n reads, 

Vk = Si I yk+m+j + ■ ■ ■ • 

^k+m+j consists of [llk+m+ji Hk+m+j+li ' ■, Hk+m+n+j)- As t/k+m+n+j = (~1) Uk+j-li it Con- 
tains {—lYyk- Thus the linear independency of solutions concludes at least 

5j;^+^) = 0, l<]<n. (24) 

We also note a trivial relation, S{1 " = (—1)"". 

In the next section, we will show that a set of solutions to the above recursion relations 
can be neatly represented by the quantum Jacobi-Trudi formula. 
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4 Quantum Jacobi-Trudi formula in Stokes multipli- 
ers 



The recursion relation (p3D can be successively solved for given initial conditions 
T^^^X), ■ ■ ■ ,T^^\\) (forward-propagation). Conversely, t^"-\X) is expressible by these so- 
lutions (back-propagation). 

To represent these, we employ notations useful in later discussions, 

T^'\X) = ^W(Ag-(^-i)/2) (25) 

Ti''\x) = i-iy+^T^'^^Xq-^^'^^/^) (= (-l)"+^52(Ag-("-i)/')). (26) 

We first state the result for the forward-propagation problem, 

T« (A) = det frf -^+^')(Ag(^+2-^-^')/2)) . (27) 

V / i<,i,j<e+i 

The rhs includes the formal extension of T{" (A), and we understand TJ^^ (X) = for 
a > n + 2 or a < and Ti°\X) = 1 . 

The proof is easy by the induction on i. The case £ = is trivial by definition. Assume 
that eq. (|25|) is valid up to i. Expanding the determinant with respect to the last column, 
we have 

tSiW = T«(Ag-^/^)T«(Agi/^)-Tf(Ag-(^-i)/2)TH(Ag) 

+ . . . + (-l)'^+iT('^)(Ag-(^+i-«)/2)T,« _^(Ag«/2) + ■ ■ ■ (28) 

Here we set T^^\X) = 0,£ < and T^^\X) = 1. 

By eqs. (^), (^) and from the assumption of the induction , the Ihs is equal to that 
in eq.(^) after A —^ Ag~^/^. Then the validity of eq.(p5D carries over into £ + 1, which 
completes the proof. 

Next consider the back-propagation problem. We regard the coupled equations (|23| ) for 
£ = 0, 1, 2, ■ ■ ■ , n, as a set of linear equations for r*^"^(A) (or T{ (A)). The Cramer's theorem 
is again applicable. Taking account of the fact that the determinant of the coefficient 
matrix is unity, we have. 



Tr(A) = det Ti^i^,+,(Ag^*+^-'^-i)/') . (29) 

V / l<,i,j<a 

It is natural to define T| (A) = 1 from this expression. 

We notice the similarity of eqs. (p7|) and (l29|) to the Jacobi-Trudi formula for Schur 



functions ||35|. It states that any complex Schur function associated to a skew Young 
diagram can be represented in terms of a determinant of a matrix whose elements are 
given by elementary Schur functions. 

We find a quite parallel result in the present problem. Thanks to the argument around 
eqs.(^) and (^), only 5"^ {(A) is of our interest. This turns out to have a compact ex- 



pression in a similar manner to the Jacobi-Trudi formula, which we call the quantum 
Jacobi-Trudi formula. 



To state this, we prepare some notations. By /i, we mean a Young diagram (/ii, /i2, ■ ■ ■ ) 
and by /i', its transpose. Consider a skew Young table of the shape n/rj = (/xi — rji, fi2 — 
?72, ■ ■ ■ ) such that fii > rji. We define a quantity associated to fi/r], 

T,/,{\) := deti<,,,<,, (T(;L,,.,+,(Ag-('^i-'^^+'^^+'"=-^--^+^)/^)). (30) 

It contains eq.(^) as a special case /i = (l,l,---,l) and f] = (p- 

A generalization of eq.(p7D also exists. To show this, consider two matrices H and E 
with entries 






Clearly, detif = detE' = 1 as they are lower triangular with all diagonal elements unity. 
The expansion of eq. (pTl) with respect to the 1st column leads to 



= -T«(A)+T«(Ag^/2)T«(Ag-V2)_j.(2)(^^(.-i)/2)yW(^^-i) 
+ . . . + (-iy+^Tt\Xq^'^'-''^/')T^l\_,{\q-'^/') + ■■■ 
which is equivalent to 

k 

Hence E and H are inverse each other. 

These two facts are sufficient for the second representation of T^/,, |^ ^ ^ , 

T,/,(A) = deti<,,,<,,(T,^^^-''^-^'^'\Ag-(^'^---^i-''^+^-+'^-i)/^)). (31) 

Eqs. (|30|) and (0) coincide with the Jacobi-Trudi formula by dropping the A dependencies. 
We have the following statement on nontrivial Stokes multipliers. 



Theorem 3 The solution to the recursion relation ( |7q j for j = 1 is given by 

Si'!iW = i-ir^Xiq^'^'-'^/'X), (32) 

where fi is a Young diagram of the hook shape fi ={£,!,■ ■■ ,1) with height k. 

The proof is as follows. In terms of T, we need to show 

T,, ( A) + T,„ ( A) = Ti'^ {W^)tP ( Ag-^/2) , 

where /^'(yu") is also the hook Young diagram //' = (£ + 1, 1, ■ ■ ■ , 1), (/x" = (£, 1, ■ ■ ■ , 1)) 
with height k (A;+l). Then the equality is verified by expanding the determinant associated 
to T^" (A) with respect to the first row. 
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It remains to show the consistency of eq. (p^ ) with eqs. ( |T8| ) {k = n + 1 and j = 1) 
and dH). 

First, consider k = n + 1 and j = 1 of eq.(|l8]). By the recursion argument in the 
above, the formal extension S'„4_2 iW^ defined by the rhs of eq. (|32D , must be zero. This is 
directly seen from the second expression (^), as T| (A) = 0, {a > n + 2). 



Next we deal with eq.(|T5|). To be precise we will show 



5f;+"+^)(A) = T«„+i(Ag(-+")/2)5,, (= i-ir5,,. (33) 



With slight re-definition of index j and the application of the relation (|21|), this is 
converted into an equivalent form, 

^(^■^^^(A) = T«„+i(Ag(-+'^)/^)5,+,v.+i (0 < J < n). (34) 

We shall divide argument into two cases, j ^ and j = 0. 

For j = 0, we apply eq.(|3^) to 5*1™ (A). Then we find from the formula (pO|) that the all 
elements of the first row are zero unless i = n + 1. This comes from eq . (plf) , Tm+j'(A) = 
for 1 < j' < n. For i = n + 1. We expand the determinant with respect to the last 
column. From the normalization in eq . (p2D and Tq (A) = 1 , it follows >S'„'^^ ^^ (A) = 
Ti||„+i(Ag('^+")/2)^ ^hich is (-1)" by definition. 

A remark is in order. The second form (^ ) leads to Sl-^^ (X) = 

(-l)"Ti!^(Ag('"-i)/2)T^("+i)(Ag3m+n/2))_ r^^ieu 

T^^\\q-^-^'^")Tt^'\Xq-") = (-l)"T«„+i(A) = 1 (35) 

must hold. Neither Tm {X) or T^^" ^(A) has poles. Thus eq.(^5D asserts that they are 
also nonzero everywhere. Consequently they are constant. We drop their A dependen- 
cies for the time being. By taking the determinants of both sides of (matrix form of) 
eq.(|3D, one obtains (T^^'^^^^ )'"+"+! = {T^l^+iT^^- Being combined with eq.(||), this 



' m+n+l) 

yields, (T{" )"* = {Tm Y'''^^- We conclude both T^ and Tm are root of unity. The 
former can be derived directly once if we assume (^ and use (0). One only has to evaluate 
Wk using asymptotic forms (|^) and (|^) in Sk+(^n+i)/2^Sk+(n-i)/2- The simple manipulation 
yields r("+'^ = (-1)", thus Ti^"^') = 1. 
For j 7^ 0, the recursion ([TsD yields, 

^(-+j+i)(^) = r«(Ar+^>ir'^(A) + 5j7+/)(A) 

= (-l)*T(*)(Ag'"+^^+(-'-i)/2)T^i|^.(Ag(™+^--i)/2) + 5J;"5'Ha) (36) 

In the present case, the first term in eq.(P3) is vanishing due to eq.(E3p. Thus 



Sl^^^^'\X) = S^^i\X). (37) 



We then try to prove eq.(|3^ by induction on the upper index, starting from the re- 
sult for Sl"f^ (A) as the "initial condition". In each induction step, however, i = n + 1 
component is indeterminate from eq.(p7D. We then use eq.(|3TD and find Sl^lii{X) = 

{—lY^^T^_\{\q^^^'')T{' {\q^ ^^^) where u!{i) and uj'{i) are some shifts irrelevant in the 
present argument. This supplies the missing pieces. As £ > m + 2 for the case under 
consideration, these are null. We then immediately check that eq.(|3^ holds. 

Thereby, we prove that nontrivial Stokes multipliers are explicitly given via formula 
(^). The fundamental quantities r'-'^-* s are not all independent. Constraints (p4D impose 
complex algebraic equations among them. 

We note useful relations among Ts, the T— system. Let T^ (A) = T^(A) for fi being a 
rectangle of height a and width i. By applying the Pliicker relation to eg. (|30|) or eq. (|3lD , 
one finds, 

Tt\Xq'/')Tt\\q-'/') = T^l\i\)Tt\i\) + Tt'''\\)Tt-'\\) , a = 1, ■ ■ ■ ,n. (38) 

The boundary conditions are T^'J^{X) = T^'|i<,<^+JA) = Ti''^'\X) = Ti^>''^'\X) = 
and T^ (A) = lifa = Oor£ = 0. The last relation comes from the null dimensionality of 
the matrix in eqs.(PIl|) and (PT]). The initial condition T^Jf^{X) = leads T^"q(A) = 0, (1 < 
a < n). Similarly we have Zm+^(A) = (1 < a, j < n). Thus the T— system constitutes 

finite number of relations with finite number of T^ (A). 

Remark that the case with ra = 1 , m = 3 and with multi-parameters has been found 
in [^ by direct calculation of recursion relations. The one parameter result has been later 
extended to arbitrary m by similar argument employed here 0. 

Summarizing this section, we have found determinant representations of Stokes multi- 
pliers. We conveniently introduce a subset of huge hierarchy, T^ (A) including (the most 
fundamental) Stokes multipliers r*^"^(A). Then the functional relations exist among them 
which, in turns, impose some restrictions on r'^'*^(A). 

We will discuss above results in view of solvable models in the next section. Especially 
we will give some follow-up of the last sentence in the previous paragraph. 

5 Functional relations in solvable models and Ther- 
modynamic Bethe Ansatz 

Hereafter we assume n odd and m = M{n + 1). It is then convenient to rotate the A for 
the differential equation (|l]) by |. We will use the same notation for the resultant A. 

In the first part of this section, we remind of some concrete results in solvable models. 
The second part is devoted to a rather speculative discussion of the possibility of the 
application of Thermodynamic Bethe ansatz to the evaluation of Stokes multipliers. 

The commuting transfer matrices play a fundamental role in studies of solvable lattice 
models and field theories [|ll|, [12], |l3l . The members in a commuting family share the same 
physical space ( quantum space) and are parameterized by the (multiplicative) spectral 
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parameter A. They are labeled by the auxiliary space of which trace must be taken. These 
auxiliary spaces are identified with irreducible modules of Yangian or quantum affine Lie 
algebra [^ Q. For Uq{An ), there exists a irreducible module Wm (A), (a < n) which 
is isomorphic to wiVa^ as a classical module. Naturally, we associate a Young diagram of 
a X m to this module, and write the corresponding transfer matrix 7^" (A). Since they are 
commutative, we consider transfer matrices on common eigenstates. Thus we sometimes 
do not distinguish operators from their eigenvalues. 

In the language of solvable lattice models, 7m" (A) is a transfer matrix of a model 
obtained by the fusion procedure. Starting from a "fundamental" model acting on Vai x 
Vai , we can recursively derive fusion models. Their auxiliary spaces are constructed from 
Vai X ■ ■ ■ X Vai by applying appropriate projectors. We utilize singular points of R— matrix 
for construction of these projectors. 

The fusion procedure is not restricted to the rectangle type. Generally, there exists 
a module iy^/r,(A), or a fusion model, parameterized by a skew Young diagram fi/rj. We 
write the corresponding transfer matrix as T^jj^^X). 



The short exact sequence of irreducible modules examined in |36, 28, 29 leads to the 



relation which we have seen in eq.(pO|) 



%/,{\) = deti<,,.<,,(T/^^^''^-"^'+'\Ag-(^i-^-^^<+^-+^-i)/2)) (39) 

= deti<,,,<,;(T^l^^_^.^,(Ag-(^'^-^^+^^+''-^--'=+^)/^)) (40) 

whereri^io(A) = r/'^<°^(A)=0. 

This assures the T— system (^) among T [^. The quantum Jacobi Trudi formula plays 



a role in several problems in solvable models [|39| , ^ |T 

The conditions (^4]) also hold. When q = exp(i ^^^ j^ ), the truncation of the space hap- 
pens due to quantum group symmetry. In view of solvable lattice models, it corresponds 
to the situation that some projectors are vanishing and fusion paths in local variables are 



lost. Consequently, no local variables can be adjacent and transfer matrices vanish p8[], 

^it (A) = T^U^) = ... = Tit (A) = 0. (41) 

The normalizations of Tq {X),T^ (A) and 7^ (A) depend on choice of quantum 
space. We have not yet found the description of the quantum space yielding Tq{X) = 
7^ = 7^ =1, except for n = 1. We assume the existence of it for general n. 

Then, on the corresponding space, the fusion transfer matrices share same functional 
relations with Stokes multipliers. 

Below we use the same symbol T^ (A) for these two cases. 

Thanks to the (quantum group) reduction,the T— system closes within a finite set of 
unknowns. 

The solution to functional relations, however, is not unique. Additional information 
on the analyticity of T^ (A) in the complex A are required for the uniqueness. 
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We start from the lattice model. Then we take the scaling limit (or the field theoretical 
limit) which has been discussed in several literatures. In the present context, we refer to 
P^ . Note that our transfer matrices are not mere sums of products of Boltzmann weights 
but are imposed overall renormalizations due to 7^ (A) = 7^ = 7^ = 1. 

Let us define an "additive" spectral parameter uhj \ = exp( . _^^^j^J . We also intro- 



duce 



Then eq.(p^) reads 

Yn '{u + i)Yf '[u-i) = , ^^, -, — -T , (42) 

where (F/ \u)Y^ = Ym (u) = by definition. 

Remember that that Boltzmann weights are regular functions of u. Then, apart from 
some renormalization factors mentioned above, T^ (A(n)) must not be singular in the 
strip Qu G [—1, 1] but may possess zeros in general. 

Suppose that T^ (A(m)) has finitely many distinct zeros in the strip {u^l} which 
depend on eigenstates. They are assumed to be off the lines '^u = ±1. 

Then eq. (^) can be transformed into coupled nonlinear integral equations by the 
standard trickfl^, ||, |T2|, |T3|, ||, gSf . 



The both sides of eq. (^21) are analytic with finitely many distinct zeros and have 
constant asymptotics. Then, applying Cauchy's theorem, we have 



\ogY}''\u) = log Zi^\u) 



+ > > / Kiu - M ) log . , ^, -, — -T du + C) 1 43) 

where C)" is a " integral" constant fixed by comparing asymptotic values of both sides. 



Zf'^u) signifies. 



z^:'\u 



zfl{u)Z^t^,{u) 



Zt\u) = n(tanh^(«-t.£)))-'^xntanh^(«-4:i). (44) 

j k 

and Zf\u) = Z^;^^\u) = zl^\u) = 1. 

Here {w) "• )} i^ the joint set of zeros (e^j = 1) and singularities (e^j = —1) of Tf" {u) 
in Q=M G [—1, 1] due to the above normalization. They stem from common factors in 
Boltzmann weights, thus they are zeros or singularities of order N [N = system size). 
Here we label them disregarding of their multiplicities. 

12 



The kernel function is easily written in the Fourier transformed form. We define the 
Fourier transformation f[k] of a function f{u) by 

-| /"OO /"OO 

^^ J -oo J —oo 

then 

^m = — L_. 

^ ^ 2coshfc 
Eq. (|43|) might be represented by the following form, in analogy with thermodynamic 
Bethe ansatz equation in CFT, 

"^ I'OO 

log r/"^) (u) = J2 ^''"(^ - "') log ^r"^ ("')^"' 

a. m /"OO 

+ E E / <i(^ - «') iog(i + (>;^'H«'))-^)c^^' + 1^1" 



(45) 



where Di' is also some "integral" constant and Y^ (m) = Y^ ('u)/Z)" (-u). 
The Fourier transformations of the kernel functions read, 

K';[k] = A^'^^m^Ak] 

-T-irrj smh.{inm{i, r)k) sinh((m — max(£, r))A;) 

siiah(mk) sinh(A;) 

2 cosh k 

and la^b = 1 if they are on adjacent nodes of Dynkin diagram for An and la^i, = otherwise. 
At this stage, one performs the field theoretical limit , system size —>■ cxd, lattice spacing 
— > 0, and sending elliptic no me —>■ under some fine-tuning condition. Precisely speaking, 
this condition depends on the regime of the lattice model. We shall skip that detail here. 
We refer H^] for detail in the case oi A{ . 



Then it can be shown that {w^" } accumulate in oo and reproduce "momentum" term 
in eq.(|5D from the first term |I], |l|, |13|, |15H|] , 



' TT "'' POO 

log f/"^ (u) = Mf^ exp(-M) + Y^ A'^'{u - v!) log zf^ {v!)dv! 

+ E E / <^(« - ^') i°g(i + {y^'\^)r')d^ + ^?^- (46) 

Z'^^\v!) consists only of {uf^} in (H. 
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The rapidity 6* in |I|] is denoted by ^m here. 
The resultant coupled integral equations can be recursively solved and yield a unique set 
of solutions to F/ (u), and then to T^ (A(m)) for given {u^l} and MJ . Thus one must 

adopt appropriate choice of these parameters to reproduce proper T^ (A(u)) including 
T^"'\\). Conversely, if two sets of functions satisfy the same relation (|38D and share the 
same {u^k} and M^""' , then they are identical in the strip. 

For solvable lattice models, the case studies on Uq{A\ \ 3) indicate that Tf''{\{u)) does 
not possess zeros in the strip '^u G [—1,1] when acting on the largest eigenvalue sector. 
That is, {ufl} = 0, a = 1, ■ ■ ■ ra for all p < m. 

This seems to be also the case with Stokes multipliers. Actually this is the case for 
n = 1. That is, the proper Stokes multipliers are reproduced by {uil} = (f),a = 1,- ■ -n 

for all p < m but different choice of M^ from lattice models or field theoretic model. 
Let me just present some argument for n > 2 case. Assume that the solution to eq. 
(^) is simple, T|" ^ = Tm = 1. As remarked above, the former equality is a direct 
consequence of (|^) and (|1T]). From our "boundary conditions", the T— system is then 
invariant under the simultaneous transformations {a,i) — > {n + l — a,m — i). This implies 

'^e W — '^m-i (-^)- By ±9f' , we mean the imaginary part of zero of T^ (A). Then a 
sum rule holds, 

f^e + Om-e -m + n + 1. 

The simplest solution, 9^'^ = a + i, is actually correct for n = a = I. Postulating this 
solution, we conclude {uy'l} = 0, a = 1, ■ ■ ■ n for all p < m also for Stokes multipliers. 

We do not expect ill-behavior for small angle sector of A; r^"'^ (A) might not be vanishing 
in the sector, or otherwise it harms linear dependency of FSS. The above choice of O^^ is 
consistent with this expectation. 

Admitting these assumptions, we have a conjecture 

Conjecture 1 The Stokes multipliers of eq.^ are given by "vacuum expectation values" 
of fusion transfer matrices related to Uq{An ). Thus they are evaluated from eq.^j^) by 
forgetting the second term in the rhs . 

The parameters M^ , Di" must be tuned correctly as in |^ for n = 1. 

To be precise, overall normalization of A for the identification is not fixed by func- 
tional relations only. For n = 1 the determination of this factor is quite straightforward 
in, H, due to the fact that only one subdominant solution exists in each Sj. The spectral 
determinant is then identified with the fusion Stokes multiplier which connects the sub- 
dominant solution on the negative real axis to the dominant solution on the positive real 
axis. Then standard WKB arguments yield the identification of parameters, especially 
M^ and the normalization of A. 

This can not be generalized straightforwardly for n > 2. Let me just comment on two 
fundamental problems. First, the meaning of eigenvalue problem is not necessary clear for 
higher order differential equations. The characterization of the eigenspace (it is the Hilbert 
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space for n = 1) is not obvious. Second, more technically, there are several subdominant 
solutions in each sector. We hope to clarify these issues in future publications. 
Finally we comment on the (1,1) component of eq.(B]). 



-(1) 



(2), 



Vk = U 'Vk+i + r^ 'yk+2 + 



■■ + Tk VkA 



■n+1- 



(47) 



The equation coincides with the Baxter's T — Q relation for n = 1 by putting a; = 
and writing yk{x = 0) = Q{Xq^^'^°'^^^'). For general n, this coincides with the "spectral 
curve" equation argued in [|^ or the characteristic equation in the quantum separation 



of variables 47 . 



6 discussion 



We have seen that Stokes multipliers associated to a n + 1 th differential equation satisfy 
the same recursion relations with fusion transfer matrices related to Uq{An ). Under the 
assumption on analyticity of some functions in the strip, we conjecture thermodynamic 
Bethe ansatz type equations eq. ( ^6]) which determine Stokes multipliers. The assumption 
, however, needs extensive numerical check, which we hope to report in the near future. 
The deformation parameter q = exp(i ^_^^ -^ ), which arises naturally in the present 
context, has a concrete meaning in the solvable models. The parameter m specifies the 



level of dominant integral weight of the local variables in the lattice model|^8|. Then 
the denominator of the exponent of q is of the form, level + dual Coxeter number. This 
combination also appears as the deformation parameter of a solvable model based on 



other affine Lie algebras |^9[. We thus expect a possible extension of the present study to 
other type of potential terms related to other affine Lie algebras. 
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Appendix 

We present the argument leading to eq. (^ . 

In terms of the vector u := {y, dy, ■ ■ ■ , d'^y), eq.(0) reads. 



d 

dx 



u = A{x)u 



/ 0, 1, 0, ■■ 


■' ^\ 


0, 0, 1, ■■ 


-, 


0, 0, 0, •• 


•, 1 


\P{x), 0, 0, ■■ 


-, oj 
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We define a new variable ^ by ^"^^ = x. Put u = MTw such that M = diag(l, ^™, ^^™, ■ ■ ■ ) 
and that T diagonalizes a matrix of the form A{x) in the above with P{x) = 1. Then the 
equivalent matrix equation is written in the form 






'w. 



and 5o = diag(exp(^|^),exp(^^^), ■ ■ ■ , 1). 
Following p^, |33], we assume w in the form, 



w 



exp(/ r^"+"7(r/)rfr7), 



P2{0 
\Pn+liOj 

and Pi(0 = 1- 

Then the substitution of the above form into eq. (^Hf ) yields, 



kPki 



(4J 



(49) 



where Pjk denotes the {i,i) th component of -8(0 and j > 2. Note that /3jfc = 0{1/^) if 

Immediately seen, Pj{i) = J2n=iPj,n^~^ is a formal solution; one can recursively 
determine pj^N by (^9|) . 

We define ho = Pn+i,n+i{oo) — /?i,i(oo), and the sector S in the ^— plane by 

3 

|arg/io + (m + n + l)argO < -vr, |0 > 17 

for a fixed positive fi . By Vr we mean the domain in (ai, ■ ■ ■ , am) space such that 

|ai|^ H la^p < r 

for some positive r. Then we naturally generalize the result in |^, |3^ . 

Conjecture 2 one can construct Pj^riO so that PiriO is holomorphic with respect to 
(^, ai, ■ ■ ■ , am) in S X T>r and d^pJ^riOi (^ — 0? l? ' ■ ■ ) admits the uniformly asymptotic 
expansions, 

for (oi, ■ ■ ■ , ttm) G I^r (IS C, tends to infinity in S. 
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Let us present some argument supporting this (by no means a proof). We put pj 



Qj + Pj,r{0 ^^d inserting this into (^Of) 



- 5^/3i.]v(0]v(0 - r'"-'^^Pi;(o (j > 2), 



«>2 



c^r 



Cj,fc = /5i,fc(l - 5j^k) + (^i./tK/^ij - /?l,l) - X] A.iPTrlO} 



j>2 



- [3l,k£r{i), U,k>2). 



The formal solution is 

qM) = / ^'"■'"(^.,1 + hkQk + 5^ A.g.g, ) X exp( ^° rr+"+^ - r7'"+"+^))rfr^. 

(50) 

where we set Cj^k = ^o "i,A: + 0j,fc and /iq = (n + l)(exp(27ri "^j|^~-^ ) — exp(27rz^^)). Now 
the coefficients in the bracket of the integrand are less than order |. 

We are interested in eq.(HD. For this, it may be enough to treat j = n + 1. Assume 



there exists appropriate choice of M, then Lemma 3 of [Q also applies in this case. The 
integral in the rhs of (|50D can be thus bounded from above , \qn+i\ < C^(|^|~^) iii a certain 
domain provided that 

|arg/iQ + {m + n + l)arg^| < — — e (51) 

where e is positive and small. By definition, Hq = ho, thus in S, g„+i tends to zero 
as 1^1 -^ oo. Hence the solution converges to the asymptotic expansion form. It may be 
shown that the limit does not depend on r as in [^, ^. The condition (|5lD coincides, 
when written in terms of x, with (H). 
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